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Abstract 

We derive and study a theoretical description for single file diffusion, i.e., diffusion in a one 
dimensional lattice of particles with hard core interaction. It is well known that for this system 
a tagged particle has anomalous diffusion for long times. The novelty of the present approach is 
that it allows for the derivation of correlations between a tagged particle and other particles of 
the system at a given distance with empty sites in between. The behavior of the correlation gives 
deeper insights into the processes involved. Numerical integration of differential equations are in 
good agreement with Monte Carlo simulations. 
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I. INTRODUCTION 

Single file diffusion, or dynamics, is the diffusion of identical particles in a one dimen- 
sional lattice, with the following features. Double occupancy is not allowed in each site of 
the lattice. In general, hard core interaction between particles is assumed. Particles can 
not jump one over the other, so that the order is preserved. The most salient feature of 
this model is that, for large times, a ta^ed particle has anomalous diffusion, with mean 
square displacement (MSD) A^x ~ t^/"^ [1]. This simple model, or different generalizations, 
has been useful to describe a wide variety of microscopic processes in physics, chemistry or 
biology (see, e. g., |2-5|). 

From a theoretical point of view, the model has motivated the development of many 
interesting results. A few examples can be mentioned. Considering a lattice spacing a and a 
homogeneous concentration per site c, Richards demonstrated that, for large times, the 
asymptotic behavior of the MSD of a tagged particle is 



c \ vr 



Is 



where V is the transition rate of a partide to an empty neighboring site. In ^ van Beijeren 

et al investigated the MSD for all times and for finite chains. Rodenbeck et al [7| studied the 
density defined in the space of all particle positions, with boundary conditions that take care 
of the hardcore interaction and, using the reflection principle, obtained an exact expression 
for the tagged particle density that slightly differs from a Gaussian distribution. KoUmann |8|] 
generalized previous results to interaction potentials, that exclude mutual passage, different 
from hardcore interaction. Lizana et al |9| obtained exact results for the density of the 
tagged particle in a finite system of length L with reflecting boundaries for three different 
time regimes: normal diffusion for short times, anomalous diffusion for intermediate times. 



and constant MSD for times longer than the equilibrium time. Manzi et al 



10| analyzed 



the role of different initial particle configurations o: 



/c-mers in the transition from normal to 



anomalous diffusion. Centres and Bustingorry [HI] described the single-file diffusion using 
a discrete Edward- Wilkinson equation. In order to have a better understanding of the 
processes that take place at atomic scale, there has been experiments carried out by Lutz 
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et al 



12| and Delfau et al [l3j. The former studied the diffusion of particles constrained 



in a circular optical trap; the latter analyzed the dynamics of particles confined in a linear 
channel of finite length. Both experiments witnessed the anomalous behavior for large times. 



More recently, Carvalho et al 



14l | investigated the effect of a repulsive Yukawa potential in 



a quasi-one-dimensional system. 

As far as we know, no diffusion equation for the density of a tagged particle has been pro- 
posed, probably because there is an important difficulty. It is necessary to take into account 
the correlation between tagged and untagged particles, since the presence of the untagged 
particles modifies the diffusion making it anomalous. This correlation, in turn, is related 
to correlations among three particles, giving rise to an infinite hierarchy of equations. An 
alternative to describe the tagged particle density with a diffusion equation is the fractional 
dynamics approach |15|. However, the fractional diffusion equation can not represent the 
two different regimes (or three for finite systems), for short and long times, that are present 
in single file diffusion. 

The novelty of this paper is the derivation of a diffusion equation for the correlation Uij 
of having the tagged particle in site i of the one-dimensional lattice and a particle (tagged 
or untagged) at distance j, with empty sites in between. In the derivation it is necessary 
to make use of some appropriate approximations in order to avoid the infinite hierarchy. 
We show that, despite the needed approximations, the equation correctly describes the two 
regimes of normal and anomalous diffusion for short and long times respectively. 



II. DIFFUSION EQUATION AND CORRELATIONS 

In the following, we will focus on a one-dimensional infinite discrete lattice. The lattice 
spacing is a, so that the position of a given site i can be obtained by x = az, Vz G Z. The 
occupation number of each site is rrij, equal to 1 or depending on the site being occupied 
or empty respectively. We define rij as the occupation number corresponding to only a 
tagged particle, while rrij corresponds to any particle. A particle jumps to a neighboring site 
with a transition rate F. Due to a hard core interaction, the jump is performed only if the 
destination site is empty. In the following, we use sans serif font for variables of a particular 
configuration of the system, and italic font for configuration averages, i.e., Ui = (rij) and 
rrii = (rrii). 
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Let us consider the tagged particle current Jj between sites i and i + 1. We can evaluate 
Jj by considering the particle configuration and the possible processes: 

M = n,(l - mi+i)r - ni+i(l - mi)T. (2) 

The first term corresponds to a jump of the tagged particle from site i to i + 1, given that 
site i + 1 is empty, and the second term represents a jump from i + 1 to i. This equation 
can be rewritten as 

J, = -(ni+i-n,)r + qir, (3) 

where = rrijnj+i — njirij+i. The densities at position x are m{x) = irti/a and n{x) = rti/a. 
Similarly, we define J(x) = (J,) and q{x) = /a?. From Eq. ([3]) we get 

J-D^^^D,, (4) 

where D = a^T is the diffusion coefficient, and we have approximated '^(^+"-y~'^(^^ ^ xhe 
diffusion equation is 

dn dJ ^ d f dn \ 



dt dx dx \dx 
The anomalous diffusion at large times can only be produced by the correlations in q. As a 

first approximation, we can decorrelate the products in q, i.e.. 



so we get 



and 



(mjiii+i) ~ (mi)(ni+i), (6) 



dn dm 

q ~ am— an——. (7) 

ox ox 



On ( On Om \ 

— - = D-— I (1 — am)— — h an-— ] (short times). (8) 
Ot Ox \ Ox Ox J 

For simplicity, we assume that the particles are uniformly distributed, so that ^ = 0. 

The tagged particle density, n, has a normal diffusion with an effective diffusion coefficient 

D' = D{1 — am). This result holds for short times when starting from random initial 

distributions, i.e. while the decorrelation ([6]) is a valid approximation. 

The crossover time between normal and anomalous behavior, tc, can be obtained by 

comparing the MSD in both regimes. For normal diffusion we have A^x = 2D't, and for 

anomalous diffusion we have Eq. ([T]), with c = am, so 



In the rest of the paper we focus on a discrete description of the system. The discrete 
version of Eq. ([5]) is 

Oti ■ 

= r [nj+i + ni_i - 2ni - {qi - . (10) 
For simphcity, we will consider a lattice spacing a = 1. 



III. TAGGED-UNTAGGED PARTICLE CORRELATION 

It is convenient to define a new correlation that is directly related to = (qj). We define 
the tagged- untagged particle correlation, the configuration average of the tagged 

particle density in site i times the particle density in site i + j (tagged or untagged), with 
empty sites in between: 



(ni(m,+i - l)...m,;+j) (11) 



with \j\ > 1. We assume average initial distributions of tagged and untagged particles that 
are symmetric respect to i = 0. In particular, in Sect. IIVI we discuss an initial random 
distribution of untagged particles, with a homogeneous density, and one tagged particle in 
i = 0. On average, the symmetry i —> —i is preserved, so that Uij = U-i^-j. We make use of 
this symmetry in order to consider only positive values of j. The relation with correlation 
Qi is 

qi = Mj+i -1 - = - (12) 

The tagged particle density is related to Uij by 

CXD 

ni = '^Uij. (13) 

i=i 

The interparticle distribution, pj, gives the probability of finding a particle at distance j 
from a tagged particle, with empty sites in between. It is given by 

oo 



-oo 



Since the tagged particles are not different from the rest in their evolution, pj actually gives 
the interparticle distribution between any pair of consecutive particles. It can be shown that 
the equation for pj , for j > 1 is 

^ = 2r + - 2p, + pi{p, - p,_i)] , (15) 



and for j = 1, 

^ = 2T[p,-p^il-p,)] (16) 

The stationary solution of these equations correspond to a homogeneous random distribution 
of particles with constant density rrii = c and is given by pf = c(l — cy~^. 

Now it is possible to derive a closed equation for Uij that is useful to obtain and to 
get a common description for the normal and anomalous behaviors of single file diffusion. 
In the derivation we have to consider all possible processes in a small time interval, 6t, that 
may modify the value of Uij at time t. Schematically, 



Uij{t + 6t) = u 



f •rx o\ ( • ^ o 

1 - Pr ■ ■ ■ - Pr • • ■ 

o \ /• o r\ 



-Pr ( ■ ■ ■ - Pr 

i+j J \ i i+j 

+Pr ( ■ ■ ■ + Pr ■ ■ ■ 

i-l i+j J \i 



o \ / • o r\ 



+Pr ■ ■ • + Pr ■ ■ ■ , (17) 

where '•' represents a tagged particle, and 'o' any kind of particle. More explicitly, we have 

u,j{t + 6t) = u,j[l-2T5t 

-T6t{l - Pr(mi_i/n,mi+j)) - T6t{l - Pr(mi+j+i/nimi+j))] 

+Ui-ij+iT6t + Uij+iT5t (18) 

+Ui+ij-iT5t{l - Pr(mi/ni+imj+j)) + Uij^iT5t{l - Pr(mi+j/njmi+j_i)), 

where Pr(mj_i/njmj_|_j) is the conditional probability of having a particle in i — 1 given that 
we have a tagged particle in i and a particle in i + j, i.e. 

Pr(mi_i/nimi+j) = — (19) 

^njrrij+j; 

We approximate ( fT9l) by 

Pr(m,_i/n,m,+,) ^ Pr(m,_i/n,) = i^^^pi^ = (20) 
We need a different approximation for Pr(mi+j+i/njmj+j). In this case we use 

Pr(mi+j+i/nimi+j) + Pr(mi_i/nimi+j) ~ 2pi, (21) 
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so that 



Pr(mi+j+i/nimi+j) ~ 2pi - (22) 



this last approximation turns out to be the simplest choice to get an equation for Uij 
consistent with the equation for pj (|T5l) when we sum over i (see the appendix). 
Using ( l20l) and ( l22l) in ( IT8l) we get the equation for Uij, for j > 1, 

dui 



+2pi{Uij - Uij^i) + ■ ^ ■ ^ (23) 

rii rii+i ' 



and the boundary condition for j 

dt 



r[2u,,i(pi-l)+M,,2 + ni_i,2]. (24) 



For random initial distributions of homogeneous concentration c with one tagged particle in 



site i = 0, the initial condition for Uij is 



"ijlt=o = '^(1 ~ ^^ifi- (25) 



IV. RESULTS 



In this section we present numerical simulation results of single file diffusion and compare 
them with the theoretical description given by Eq. ( |23l) . 

We consider a one dimensional lattice of L sites with periodic boundary conditions. The 
lattice is large enough in order to avoid boundary effects during the time of the simulations. 
The initial condition is a homogeneous random distribution of particles with density c, with 
one tagged particle in site i = 0. In this case, the interparticle distribution pj does not evolve 
in time, and is pj = pf = c(l — c)^~^ for all times (note that this is the same distribution as 
for the case of non-interacting particles). We are interested in the evolution of the tagged 
particle density, its mean square displacement and the correlations. 

Fig. [1] shows the tagged particle density, n^, for different times in the anomalous regime. 
We found a good agreement between Monte Carlo simulations (dots) and numerical inte- 
gration of Eq. ( j23l) (curves). The results can be approximately fitted with a Gaussian 
distribution whose variance increases with time as t^^^. 
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FIG. 1: Tagged particle density rij, against position for different times in tlie anomalous regime: 
t = 10^ (□), t = 10^ (0) and t = 10^ (o). Concentration c = 0.7 was used. Dots correspond to 
Monte Carlo simulations in a lattice of length L = 1000; number of samples: 100000, and F = 1/2. 

n 

The curves are the results of numerical integration of Eq. ()23p related to rij by (I13p [la ] . Each of 
the about 700 particles of each sample was considered as one tagged particle in order to reduce the 
statistical error. 

In Fig. [2]we show the correlation qi for different times. Eq. (jlj) shows that this correlation 
is part of the particle current. The physical meaning of is that of a current opposite to 
the free diffusion current, proportional to the concentration gradient. The net effect is a 
reduction of the total current due to the interaction with the surrounding particles. For 
short times (normal regime) this effect is manifested in a reduction of the effective diffusion 
coefficient: D' = D{1 — c). For large times (anomalous regime) the effect of qi is more 
complex and it is not simply a reduction of the diffusion coefficient but a qualitative change 
in the temporal behavior of the MSD of rii that now increases as t^^^. 

Fig. |3] shows the MSD against time in log-log scale. For short times we have the normal 
diffusion with A^x = 2r(l — c)t. For large times, we have anomalous diffusion with A^x ~ 
t^/^ [Eq. (II])]. Results from numerical integration of Eq. (!23|) agree with Monte Carlo 
simulations. The only effect of the approximations used in the previous section is that the 
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FIG. 2: Correlation qi = ii_j„i i — Uj^i for different times. Same parameters as in Fig. [TJ 



MSD of the theoretical results is slightly smaller than the one of the Monte Carlo simulations. 



V. CONCLUSIONS 



We derived a theoretical description of single file diffusion that is different from previous 
approaches in one important aspect: it allows not only the derivation of the tagged particle 
density but also of the correlations between tagged and untagged particles. We showed that 
the correlation given by qt = u^i-i^i — Ui^i plays an important role in the diffusion equation 
for the tagged particle density, Eq. ([5]). It is equivalent to a current, opposite to Pick's law 
current, that represents the impediment to diffusion generated by the surrounding particles. 
Knowing the behavior of correlations allows a deeper understanding of single file diffusion. 

Despite some approximations needed to avoid an infinite hierarchy of equations, the equa- 
tion obtained for the correlation Uij correctly describes the normal and anomalous regimes 
for short and long times respectively. This is particularly clear for the time dependence of 
the MSD of a tagged particle, for which the behaviors A^x ~ t (short times) and A^x ~ t^/^ 
(long times) were obtained. 

Finally, let us note that, starting from a random distribution of particles, pi = pf = c 
for all times. As the number of collisions per unit time of a tagged particle is proportional 




FIG. 3: Mean square displacement against time in log-log scale. Parameters as in Fig. [TJ Dots 
correspond to Monte Carlo simulations and the curve to numerical integration of Eq. ()23p related 
to rii by (113p . Dashed lines represent the asymptotic behaviors for short [2r(l — c)t] and long [Eq. 
([T|] times, i.e., normal and anomalous regimes respectively. 

to pi, this number does not depend on time. Then, the different behaviors of the MSD at 
short (A^x ~ t) and long times (A^x ~ t^^^) are not due to the number of colhsions as 



it was aheady pointed out in [10||. On the other hand, the combination of the number of 
colhsions per unit time between a tagged particle at site i and another particle at site i + 1 
(oc Ui^i) and the number of collisions per unit time between a tagged particle at site i + 1 and 
another particle at site i (oc Ui+i-i = plays a crucial role, because the difference 

Qi = 1 — 1 appears in the evolution equation of the tagged particle density rij, Eq. 

( ITO|) . In other words, the frequency of collisions that take place at any site is irrelevant in the 
understanding of normal and anomalous behaviors, but a specific combination of frequency 
of collisions of particles in sites i and i + when we take into account in which site the 
tagged particle is, plays a central role in the behavior of the MSD. 
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Appendix 



Here we explain in more detail the approximation of Eq. (!22|) . If we take the sum over i 
in Eq. f lTS]) . we get 

E ^ = 2r [P.+i + P,-i - + (a, - a,_i)/2] , (26) 

i=— oo 

where 

oo 

= E Mij[Pr(mi_i/n,mi+j) + Pr(mi+j+i/nimi+j)] (27) 

j=— oo 

Comparing ( 126|) with Eq. ( !T5|) we obtain 

oo 

= 2pipj = E Mi,j2pi. (28) 

i=— oo 

Now, using ( 1271) and ( 128|) . we find that the simplest choice to approximate the sum 
Pr(mi_i/nimj+j) + Pr(mi+j+i/nimj+j) is the one given in (j2Tl)- 

The approximations are also consistent with the equation for rii (fTOj) when we take the 
sum over j of Uij. 
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